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Abstract With the completion of Run I of the CERN Large Hadron Collider, particle physics 
has entered a new era. The production of unprecedented numbers of heavy-flavoured hadrons 
in high energy proton-proton collisions allows detailed studies of flavour-changing processes. 
The increasingly precise measurements allow the Standard Model to be tested to a new level 
of accuracy. Rare b hadron decays provide some of the most promising approaches for such 
tests, since there are several observables which can be cleanly interpreted from a theoretical 
viewpoint. In this article, the status and prospects in this field are reviewed, with a focus on 


precision measurements and null tests. 
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1 Introduction 
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Among the most distinctive features of the Standard Model (SM) of particle 
physics is the organisation of “flavours” of quarks and leptons. Flavour-changes 
can occur only through the charged current weak interaction, so transitions be¬ 
tween fermions of the same charge can only occur through loop processes 11 . The 
probabilities of different transitions are governed by the elements of the appropri¬ 
ate fermion mixing matrices. In particular, the fact that the Cabibbo-Kobayashi- 
Maskawa (CKM) quark-mixing matrix |2|3j is found to be approximately diagonal 
suppresses generation-changing transitions. 

Consequently, processes involving flavour changes between two up-type ( u , c, 
t ) or between two down-type ( d , s , b ) quarks, i.e. involving a flavour-changing 
neutral current (FCNC), occur only at loop level and are predicted to be rare 
within the SM. Decays of b hadrons into final states containing a photon or a 
dilepton pair (e + e _ , /i + /r _ ) are of particular interest, and are the main topic 
of this review. The rates and various kinematic distributions as well as CP 
asymmetries, and other properties, of such rare decays can be predicted in the SM 
with low theoretical uncertainty, while the measured quantities may be affected 
by physics beyond the Standard Model (BSM), also referred to as “New Physics” 
(NP). Comparisons of the predictions with the measurements therefore provide 
sensitive tests for BSM contributions. 

This reason for interest in b hadron decays has been known since before the dis¬ 
covery of the b quark itself, and rare decays have been investigated by a number 
of experiments. The discovery of the b —> sy process by the CLEO experiment |4j 
has been followed by increasingly precise determinations culminating in results 
from the BaBar J5,|6| and Belle |7| experiments. The consistency of these mea- 
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surements with the latest theoretical prediction [8] provides strong constraints 

on BSM models. Among the many other important results from the B factory 
experiments, the first evidence for the B + —> t + u t decay J9 - 131 is particularly 
germane to this discussion. The overall picture is one of consistency with the 
SM, but at a level of precision that mandates further experimental investigation. 

The Large Hadron Collider (LHC) at CERN [l4] provides the opportunity to 
make the next leap in precision. Its high energy proton-proton collisions give 
a large cross-section of 0(100 fib) 1151 for production of b quarks. Due to the 
high luminosity delivered by the LHC, the decay products of the b hadrons that 
emerge from fragmentation are recorded in sufficient quantities to allow studies 
of rare decays at the ATLAS [16|, CMS |17j and LHCb (l8j experiments. For 
ATLAS and CMS, which instrument the central region of pseudorapidity, the 
online selection (“trigger”) requirements select only b hadron decays that contain 
a dimuon pair. The LHCb detector, however, covers the forward region where 
b production peaks, and is designed to enable a broader range of b hadron decays, 
including those containing a photon or a dielectron pair, to be triggered; to 
achieve this LHCb must, however, operate at a lower instantaneous luminosity 
than the other experiments. In the LHC Run I data-taking period, corresponding 
to the calendar years (2011) 2012, when collisions were at centre-of-mass energies 
of (7) 8TeV, ATLAS and CMS each recorded approximately (5) 20 fb _1 , while 
LHCb collected around (1) 2fb _1 . These data samples contain unprecedented 
yields of numerous interesting rare decays of b hadrons, as will be discussed. 

The focus of this review is the impact of the results, in the field of rare decays 
of b hadrons, from Run I of the LHC. This includes discussion of relevant results 


from other experiments, and a forward look to Run II and beyond. In order to find 
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small deviations from the SM predictions, it is essential to aim for high precision, 


and therefore observables that can be both cleanly predicted and well measured 
are of greatest interest. Such observables include relative and absolute rates, 
properties of kinematic distributions, and CP asymmetries of decays involving a 
dilepton pair or a photon in the final state. Certain processes that provide null 
tests of the SM, for example lepton flavour or lepton number violating decays, are 
also relevant in this context. This selection of observables does not by far include 
all interesting measurements in quark flavour physics, or even in B physics. The 
interested reader is referred to reviews covering CP violation in hadronic b hadron 


decays 19-21 , the B ^ system 22 , D physics 23 , rare kaon decays 24 25 and 


top quark properties 26, 27 . An earlier review on rare b hadron decays can be 


found in Ref. 28 


The remainder of the review is organised as follows. In Sec. [2] the theoretical 
framework is set out, while in Sec. [3] the experimental results are summarised. 
These two aspects are brought together in Sec. [4] to enable interpretation of the 
results in the context of the SM and BSM theories. A brief summary concludes 
the review in Sec. [U 


2 Theoretical framework 

The main challenges to develop the theory of rare ^-decays in the LHC era are to 
improve precision of the predictions, and to perform and refine interpretations of 
the data in order to map the borders of the SM and possibly detail BSM features. 
The focus is on exclusive decays of b hadrons, including B J? mesons and b baryons. 
As regards the predictions, theory greatly benefits from the determination of cru¬ 
cial input such as quark mixing and masses from earlier experiments, and from 
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maturing heavy quark methods for precision calculations of decay amplitudes. 


The latter are based on the separation in energy scale between the mass of the 
6-quark and the energy scale of QCD (mj Aqcd)- It is also possible to con¬ 
struct observables that are intrinsically robust against theoretical uncertainties, 
and hence provide precise tests of the SM with clean interpretation. In this Sec¬ 


tion the framework for these tests is outlined. In Sec. 2.1 the effective low energy 
Hamiltonian, whose induced couplings (the so-called “Wilson coefficients”) can 
be used to describe the phenomenology of a wide range of decay modes, is in¬ 
troduced. The status and recent advances of methods to determine QCD effects, 


in particular in exclusive 6 —> s(. + £ decays, are briefly reviewed in Sec. 2.2 


In Sec. 2.3, the optimised observables that are investigated experimentally are 
introduced, together with a discussion of consistency checks based on symme¬ 


try relations. Finally, in Sec. 2.4 several explicit BSM theories, are discussed as 
examples of how deviations from the SM may appear in experiments. 


2.1 Model-independent analysis of b — > s transitions 

The large masses of the W±, Z and top quark compared to that of the beauty 
quark allow the construction of a low energy effective field theory for |AH| = 
\AS\ = 1 transitions, with Hamiltonian 

HeS = ~ 4 7T VtbK £ £ C ^)°M > (!) 

v i 

where Gf is the Fermi constant, V t j are CKM matrix elements and a e is the 
fine structure constant. The Ci(/i s ) are Wilson coefficients corresponding to 
local operators with different Lorentz structure, Oi(/r s ). The operators and 
their Wilson coefficients are evaluated at the renormalisation scale /x s . Dou¬ 
bly Cabibbo-suppressed contributions to the Hamiltonian oc V ub V* s have been 
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neglected. Details of the effective Hamiltonian of Eq. can be found, for ex¬ 


ample, in Refs. 29,30 


The following local operators are important for rare radiative, leptonic and 
semileptonic b hadron decays 


0 7 

= m ^ so^P R bF llu , 

0 7 

= m e b sa^P L bF lw , 

08 

= gs^sa^P R T%G“ u , 

0^ 

= fjs sa ,iv P L T°bG' 

09 

= sy^P L blYi, 

09 

= s^p R bhn, 

0,0 

= sy^Pibty^y^l , 

0io 

= s'y fl P R bi'y^'y 5 i. 




(2) 


Here, p l/r = (1 =p 75)/2 denotes a left/right handed chiral projection, T a rep¬ 
resents the generators of QCD, and (G® u ) is the electromagnetic (chromo- 
magnetic) field strength tensor. The chirality-flipped operators 0\ correspond to 
right-handed couplings and are obtained from the O- by replacing Pl «->■ Pr. The 
left-handedness of the charged current interaction means that the Wilson coeffi¬ 
cients C[ corresponding to these primed operators are suppressed by 0(m s /mb) 
in the SM. 

The Wilson coefficients Cp can be determined from measurements of observ¬ 
ables in various different b hadron decay channels. Among the operators of 
Eq. Q, radiative b hadron decays receive contributions from Op and purely 
leptonic decays from Op . Semileptonic b —> sl + l~ decays receive contributions 
from all of Op , 0,p and 0\ q . The b —> d'y and b dl + i~ transitions are treated 
in an analogous way, but are further suppressed by V tb V bd as opposed to V tb V t * in 
Eq. 0- Consequently, in b —> d transitions, CP violation effects are generically 
larger because V ub V* d is of comparable magnitude to V tb V t * d , though the stronger 
GIM-suppression of the V b V*, term limits the size of any CP asymmetry. In the 
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SM, scalar and pseudoscalar operators 

O s = sP R bU , Q' s = sP L bU , 

(3) 

O p = sP R bhd, O'p = sP L bhd 

are highly suppressed due to the small masses of the leptons, and can be safely ne¬ 
glected even for decays involving r leptons. Contributions from tensor operators, 


O t = souvbto^i, O T5 = sGnvbla^^l , 


(4) 


are also negligibly small in the SM. 

The Wilson coefficients at the weak scale are obtained from matching ampli¬ 
tudes of the full electroweak theory onto H e s- Below the W mass, the Wilson 
coefficients follow renormalisation group evolution assuming SM dynamics |29|. 


The values at = m^ are 31 


= -0.3 , C 9 S1V1 = +4.2 , = -4.2 


-(SM 


SM 


(5) 


Comparisons of the measured values with the predictions provide sensitive tests of 
the SM. BSM theories can modify the Wilson coefficients of H e s, Eq. (jTj) , includ¬ 
ing those of operators not present or suppressed in the SM, cP = Cp SM +(7^ NP . 
The number of possible new operators, at dimension six, is large and includes 
scalar, pseudoscalar and tensor operators. If the BSM physics does not couple 
universally to leptons then sets of operators need to be considered separately for 
the different lepton flavours. New operators can also, in principle, induce lepton 
flavour-violating processes that are forbidden by accidental symmetries of the 
SM. Operators can also be associated with new sources of CP violation, making 
their Wilson coefficients complex-valued. The large number of possible operators 
is intractable for a fully model-independent analysis. However, certain experi- 
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mental signatures that can only be explained by particular extensions to the SM 


operator basis allow for a simplified analysis. 

Recently it has become customary to rewrite the semileptonic operators of 
Eq. 0 in a basis with left- and right- projected le 

^ll = (^9 — ^io)/2) C’lr = 

^RL = (^9 — ®lo )/%! C’rr = 

where 

C’ll = ^9 Cjo) C'lr = 

C RL = C' 9 -C[ 0 , C RR = 

This basis change is useful in frameworks where BSM physics at a high mass scale 

respects the SU(2 )l part of the SM gauge symmetry group, resulting in a simpler 
structure. For instance, instead of fitting the two parameters C 9 and C 10 , , if it is 
assumed that BSM physics contributes to O ll only, the constraint C 9 + C 10 = 0 
can be used. In addition, SU(2)L-relations between 6-decay observables and top 
physics can be obtained [34]. 


sptons 32,33 


(O 9 + O 10 )/ 2, 
(O' 9 + O[ 0 )/ 2, 


( 6 ) 


^9 + Qo 3 

c'+cv 


(7) 


2.2 Non-hadronic b hadron decays in QCD 

Semileptonic heavy-to-light 6 hadron decays such as B° —> K*°n + fr~ have partic¬ 
ularly interesting phenomenology. These decays have sensitivity to electroweak 
physics in two kinematic regimes: at low invariant dilepton mass-squared ( q 2 ), 
where the emitted hadron is energetic (E 3> Aqcd in the 6 hadron rest frame), 
QCD factorisation (QCDF) applies [35]; at high invariant dilepton mass, the re¬ 
gion of low hadronic recoil, where q 2 = 0(m 2 ), an operator product expansion 
(OPE) in 1 /mb applies |36|. These different kinematic regimes are indicated in 
Fig-0 In both regimes, the heavy-to-light decays can be predicted systematically 
from QCD. The methods to do this are now commonly employed, and in view 
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Figure 1: Cartoon of the differential decay rate of B° —> K*°fi + fr~ as a function 
of q 2 . At very low q 2 (near maximal Ej<*), the virtual photon contribution 
associated to C^p dominates. As q 2 increases there is a region from 1 < q 2 < 
6GeV 2 /c 4 where interference between O-j and Og becomes large giving excellent 
sensitivity to NP in Cg. At intermediate q 2 , the spectrum is dominated by the 
narrow J/0 and ip(2S) resonances. At large q 2 (small Ek*) contributions from 
broad charmonium resonances, above the open charm threshold, can be treated 
with a local OPE. 

of the experimental situation control of uncertainties becomes central. The dom¬ 
inant systematic uncertainties are parametric uncertainties from the hadronic 
transition form factors, 1 /nib power corrections (at low q 2 ), and backgrounds 
from cc resonances above the open charm threshold (at high q 2 ). In view of these 
issues it is mandatory to study the low and high q 2 regions separately, and it has 
become conventional to perform analyses in finer bins of q 2 . 

The transition form factors for heavy-to-light decays can be computed using 
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the method of light cone sum rules (LCSR) if the final state hadronic system is 

energetic. Determinations for B —> K and B —> K* form factors can be found in 
Refs. 37j381 and |39|, respectively. Recently it has also been possible to determine 
the same form factors from lattice gauge theory. Lattice calculations are only 
applicable when the hadron is almost at rest in the b hadron decay frame, i.e. 
at low recoil. The LCSR and lattice results are therefore complementary to each 
other, covering different kinematic regimes. Unquenched lattice determinations 


can be found for B —> K 40 , B —> K* and B® —> <p 41 , and A? —> A 42 form 


factors. All lattice and most LCSR determinations assume that the light-particle 
being produced is both narrow and stable. This may not be a good approximation 
in some cases, in particular for the B —^ K* transition due to the large width 
of the K* resonance. There are prospects for an improved treatment in future 


lattice calculations 43 


Issues of 1 /mb corrections exist for small q 2 only because the power corrections 
at large q 2 are parametrically suppressed, bringing them to the few percent level. 
The topic of power corrections has received a great deal of recent attention |44j- 
46 . Eventually it will be possible to determine the corrections from data [47], or 


to subject them to consistency checks, as discussed in Sec. 2.3 


Amplitudes for rare semileptonic decays also receive contributions from the 
more prevalent b hadron decays to final states containing charmonia, through the 
quark-level transition b —> ccs , where the cc resonance subsequently decays into 
dileptons. In the effective theory, such contributions are induced by four-quark 
operators 


O, = £sy M P L bcr-P L c, 0 2 = fs^P L cc^P L b. 


47T - 


( 8 ) 


These operators arise predominantly from tree-level W exchange and have large, 
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order one Wilson coefficients at the 6-quark scale: C\ ~ —0.2, C 2 ~ 1.1. As 
shown in Fig. [lj the most prominent effect of charmonia in semileptonic decays 
is the narrow resonance peaks at q 2 = ; these q 2 regions have to 

be removed in experimental analyses. In principle, however, the presence of 
charmonia affects the entire q 2 region. In radiative decays and at large recoil 


all charmonia are off-shell and are suppressed 37, 48 . At high q 2 a number of 


broad cc resonances can contribute 149,50 . Such structure has been measured 


quite precisely in the B + —> K + fi + fi decay |51|, as shown in Fig. [ 2 ] Under 
a naive factorisation assumption, this structure can be compared to that from 


e 1 e 


hadrons 152,53 using dispersion relations 491. Such a comparison 54 


revealed a dramatic deviation from expectation. Although naive factorisation is 
not expected to be valid at high precision, this surprising difference needs to be 
understood in order to maximise the sensitivity to NP. 

For completeness, it should be noted that there are also contributions from four- 
quark operators, referred to as “QCD-penguins”, with flavour structure O 3..6 ~ 
sq pPisbYluds cbQ'y IJ, PL,Rq- Their SM Wilson coefficients at the scale are 
small, 0( 1CU 2 ), and therefore their effects in radiative and semileptonic b hadron 
decays are in general small. Contributions from light resonances such as the 
4> meson, mediated by these operators, will however become important as the 
precision improves. 

The OPE does not describe the resonance contributions locally in the ampli¬ 


tudes 55 . It is, however, expected to capture their effect after integrating over a 
sufficient range of q 2 . It is therefore important to investigate the optimal binning 
for precision BSM searches. As the OPE does predict universality between cer¬ 
tain transversity amplitudes, cancellation of the effects of resonances is expected 
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Figure 2: Background subtracted dimuon mass distribution of B + —> K + ^, + ^l~ 
candidates f5l] . Broad resonant contributions from the decays B + — » -0(377O)/F + 
and B + —» ^(4160 )K + are clearly visible. 


in particular ratios 


56 


2.3 Optimised observables and symmetry relations 

Semileptonic b —>• s£ + £~ decays with a vector meson in the final state provide 
a particularly rich set of observables that can be accessed through the angular 
distribution of the particles in the decay. The decay rate for these semileptonic 
processes can be expressed in terms of decay amplitudes for the vector meson 
corresponding to different transversity states, yl^’ R , 7l RR and Aq' R . The axial- 
vector coupling through Cio allows to distinguish between amplitudes with a left- 
and right-handed chirality of the dilepton system. For completeness, it should 
be noted that a fourth amplitude (At) exists corresponding to the spin-0 £ + £~ 
configuration. The effect of this amplitude is however suppressed by the small 
lepton mass. In BSM models with (pseudo)-scalar operators the effect of At can 
be enhanced, and a further amplitude, As, appears. In the presence of tensor 
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operators, six additional transversity amplitudes can enter 57 . 


It is then possible to construct angular observables that are related to these 
decay amplitudes, for example the fraction of longitudinal polarisation of the 
vector meson, 

K | 2 + K ! 2 


F l = 


(9) 


\A^ + |A*| 2 + |Ajf + \Af \ 2 + \A \| 2 + |A R | 2 ' 

One of the most widely discussed observables is the forward-backward asymmetry 
of the dilepton system, 


Afb = 


Re ( AM* - AM* 


( 10 ) 


3 |Aq | 2 + |Ag-| 2 + |A|jf + \Af | 2 + \A \| 2 + |A*| 2 ’ 
which arises from the chiral nature of the coupling to the leptons and flips sign at 
q 2 « 4 GeV 2 /c 4 


35 


58 


59 


due to interference between the photon dipole (O^) 


and vector operators (Og^). 

It is also possible to build sets of “optimised” observables from the decay 
amplitudes in which the leading form factor uncertainties cancel. The design of 
such observables is based on the 1 jm^ expansion. In both kinematic regions of 
interest (low q 2 and high q 2 ) the transversity amplitudes receive their leading 
contribution from a factorisable form factor term 


A^’ R = ft x C7 L,R + non fact., * = 0, ||, _L , 


( 11 ) 


where the fi are corresponding transversity form factors and contain long-distance 
(QCD) information only. The C L,R on the other hand denote short-distance 
coefficients, sensitive to electroweak scale physics. They are composed of Wilson 
coefficients of the semileptonic four-fermion operators as well as contributions 


from four-quark and dipole operators (see Sec. 2.1). 

At low q 2 , the form factor relations imply that /y = f± + 0(l/mi ,), such that 
A RR ~ — A I | ,R if C = 0. This relationship allows to construct SM null tests with 
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sensitivity to right handed currents, for example [601 

|^| 2 + |A5| 2 -|A]f| 2 -|Af|2 

T \A\ | 2 + |^|2 + |4|2 + |^R|2 


15 


( 12 ) 


which is expected to be very close to zero in the SM. It is also possible to build 
other clean observables from bilinear combinations of amplitudes such that fu, 


f l and /o cancel at leading order in 1/mb. The P'-family 61 is a good example 
of a set of these clean observables, e.g. 


V2 


+« i 2 ) (141 2 +i^jr 2 +lyj 2 +i^i 2 

_ Re (C 0 L Cj - C„ R C») _ 

/(|C„ L I 2 + |C 0 L P) (icjy + |C-L|= + |C-L|= + |c{f 


(13) 


At large q 2 , the OPE predicts a different kind of relationship 


4 L,R _ 


= /|j, 0 x^’ R , A^ = f ± C: 


L.R _ 


-tL,R 

'+ ’ 


(14) 


C_’ = C 


+ 


if and only if C' = 0 . 


It follows immediately that the ratio Aq’ R /aJ[’ R is short-distance-free, i.e. is inde¬ 
pendent of Ci. This feature allows to extract the form factor ratio fo/f\\ directly 
from data. If there are no right-handed currents, C' = 0, the form factor ratio 
/_l_//i| can be extracted in this kinematic region from Alp . This information 
about the hadronic system, can in turn be used to provide better control of un¬ 
certainties on other observables. It is also possible to derive form factor free 


observables at large q 2 , such as the family 


62 


At zero recoil, when q 2 = (mg — mjj.) 2 , there are exact relationships between 
the amplitudes due to ambiguity of the direction of the K* in the B rest frame. 
At this kinematic end-point, A|(’ R = 0 and A.jp R = —\/2Aq' R [561. The endpoint 
predictions can be compared to data, as shown in Table [1} to provide a consistency 
check. Parity selection allows to extend the prediction to the vicinity of the 
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Table 1: Data on B —> K*£ + £~ and B® —» (j>£ + £~ in the endpoint-bin q 2 E 
[16,19] GeV 2 (LHC experiments) or otherwise q 2 E [16 GeV 2 , <? 2 max ] versus the 
endpoint prediction. Note that S3 = 1/2(1 — Fl)A^\ Adapted from Ref. 156 . 



Fl 

s 3 


S 7 

Ps/Afb 

S 8 /S 9 

Endpoint 

1/3 

- 1/3 

\/2 

0 

x/2 

-1/2 

B K*e+£~ 

0.38 ±0.04 

-0.22 ±0.09 

o-7o ±S:» 

n 15+ 016 

u - 10 -0.15 

1.63 ±0.57 

-0.5 ±2.2 

B® <j>£ + £- 

0 16 + 0 - 18 
U - 1D -0.12 

o-i9±8:l? 

- 

- 

- 

- 


endpoint, where observables such as Afb and P$ vary linearly with the modulus 
of the K* three momentum in the B centre-of-mass frame, and the slope provides 
a test of the SM. 

There are also more general relationships between the angular observables Ji 


of the B -± K*£ + £~ decay, defined in Sec. T3 due to the composition of the 
observables in terms of pairs of transversity amplitudes. These relationships, 
which are valid over the whole range of q 2 , serve as a further consistency check 


of the experimental results 63 


2.4 Benchmarking NP 

An important motivation when building BSM models is the origin of electroweak 
symmetry breaking and the stabilisation of the weak scale, i.e. the so-called “hi¬ 
erarchy problem” (see e.g. Ref. 1641). Among the most promising approaches are 
models that invoke supersymmetry, extra dimensions, new strong interactions 
or combinations thereof. Generically, in these models flavour violation beyond 
the CKM matrix is induced. Severe constraints on masses and couplings must 
be imposed for the NP to maintain contact with the electroweak scale: with¬ 
out any flavour suppression (which could be SM-like, with GIM- and CKM-like 
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effects) then NP is pushed up to scales as high as 10 5 TeV 


65 
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. In turn, the non¬ 


observation of BSM effects in the flavour sector provides important information 
about NP. 

Minimal flavour VIOLATION Standard Model extensions can be classified 
according to their amount of flavour violation. The term minimal flavour vi¬ 
olation (MFV) is used for models where flavour is broken in an SM-like way. 
Formally, within MFV the spurion fields which break the flavour symmetry that 
would be present in the SM in the absence of quark masses correspond to the 
SM Yukawa couplings [66 . Hence the flavour violation can be parametrised in 
terms of quark masses and CKM elements, which are known parameters of the 
SM. Any deviation from MFV corresponds to NP. The MFV paradigm provides 
an attractive way to resolve the tension between the expectation that the NP 
scale should be 0(1 TeV) due to naturalness arguments, while limits from FCNC 
processes assuming generic NP couplings point to a much higher scale. Still, 
viable non-MFV models exist with BSM around the TeV scale or higher. 

The MFV framework can be tested through CP violating observables as well 


as rare decays. Since, as discussed in Sec. 2.1, the effective Hamiltonians for 
b —> d and b —> s transitions share the same structure, ratios of b —> d and 
b —> s processes provide powerful tests of MFV. Generically, such ratios can 
be predicted in the SM to be equal to \Vtd/Vt s \ 2 ~ 1/25 modified by hadronic 
matrix elements and phase space factors, while in non-MFV models they can 
take very different values. One important example is the ratio of B° and B [? 

and B(B + 


dimuon decay rates 67 , but B(B [> —> py)/B(B u —> K*° 7) 


68 


69 


7 t + C + £ )/B(B + —> K + C + £ ) |70| exhibit similar features. To reach high preci¬ 
sion, good control of SU(3)-breaking in the hadronic matrix elements is required. 
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The ratio B(£°-> K* 0 £+e~)/B(B° -> K*°£+£~) may provide a complementary 


approach to control such uncertainties. 


Model building and SIMPLIFIED MODELS In addition to modifying FCNCs, 
NP models typically predict new particles that can be searched for at the LHC. 
As evidence for these new particles has not been found in Run I, limits on the 
masses of these particles have been pushed into the TeV range j71 . To keep the 
relation to the weak scale, without invoking any fine tuning of parameters, new 


directions of model building have appeared, e.g. Ref. 72-74 


Instead of building models that are complete up to the GUT or Planck scale, 
it is common to consider simplified models. These usually comprise the SM plus 
one new sector with a rather small number of new parameters, making them 
predictive and easy to constrain. The choice of new sector that is added is either 
made by theoretical prejudice or is driven by a need to explain a deficiency of 
the SM or a discrepancy between SM predictions and data. Two such simplified 
models, motivated by current hints of discrepancies in the B —> K^£ + £~ data, 
are Z-penguins and leptoquarks. 

A Z-penguin is a FCNC involving a neutral external field that originates from 
a U(l) gauge interaction. In the case of the SM, Z-penguins arise at loop-level 
and are induced by the weak interaction. Modifications to the effective couplings 
arise generically in most SM extensions [75 . For the b —>• s transition, 


C z = Z^(g^ b s'y fi P L b + gf b s'y fi P R b) + h.c., (15) 

where the couplings gf b R , which can be related to the Wilson coefficients Cg \ 0 , 
are generically complex. If the couplings to the leptons are SM-like, the vec¬ 
tor current coupling is suppressed relative to the axial-vector one by a factor 
14 sin 2 9w — 11, where 9\y is the weak mixing angle, and the main contribution is 
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through c[ q. With U(l)-extensions of the SM ( Z' models), however, the contri¬ 
bution to dp can also be sizable. An example of a Z' model is a gauge extension 


to t — fi lepton number 76, 77 . A survey of the parameter space of Z' models 


can be found in Ref. [78]. 

Leptoquarks are bosonic particles that carry one lepton and one quark flavour 
quantum number. They can be spin one but are more commonly assumed to be 
scalar particles (4>) which have Yukawa-like couplings (Xp, A^) to the (left- or 
right-handed, respectively) SM fermions, 


£ L Q = -Xp <t> (qP L l) - A$ ct> (, qP R Z ) + h.c. 


qt { 


(16) 


Tree-level ^-exchange induces processes such as b —> ( s , d)ZZ, which, depending 
on the handedness of the interaction, results in a modification of the semileptonic 

Leptoquarks can also provide a natural explana- 


Wilson coefficients Cg 10 


33 


tion for non-universal couplings to leptons in b —> sZ + Z~ processes. Generally, 
leptoquarks also induce lepton flavour violation, requiring an extension to the SM 
operator basis. Limits on B{fi ± -» e±j) and B(B^ —> strongly constrain 

the couplings involving electrons and muons, but the parameter space for other 
couplings remains viable. With benchmark masses between 1-50 TeV, motivated 


by the hint of lepton nonuniversality discussed in Sec. 3.4, it would be challenging 
to observe directly produced leptoquarks at the LHC but effects could be visible 
in rare decay processes such as B —» /\r ± ^ =F and 


79 


81 
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3 Status and prospects of measurements 
3.1 Dilepton decays 


The leptonic B®—> and B° —> n + decays are exceedingly rare in the SM. 

In addition to being loop- and CKM- suppressed, the decay of a pseudoscalar B 
meson into a pair of muons has significant helicity suppression in the SM. The SM 


values of the time-integrated branching fractions 82 can be expressed as 83,84 




) = 


\v; b V tq \ 2 Gl o? e M B 0 Mlfl, 

(s) ^ 


167T 3 r g H 


1 


4 M 2 

1 “ M|f l Clo(m ‘)| 2 + 

B {s) 

(17) 


where the ellipses denote subleading terms. In the above, M denotes the mass of 
the particle in subscript and T g n is the total width of the heavier of the two mass 
eigenstates in the HjU-HjU system. While those quantities are all well-known 
from experiments, the decay constant f B o must be determined from lattice QCD. 


(*) 


Using values of f B o obtained by averaging different calculations 


>) 


85 


88 


and in¬ 


cluding higher-order QCD and EW corrections 189,90 , the latest SM predictions 
are [84| 


B{B° S -> //+//“) = (3.65 ± 0.23) X 10" 9 and 

(18) 

B{B °->- n + ir~) = (1.06 ± 0.09) X 10“ 10 . 

The uncertainties on the SM predictions mainly come from the knowledge of the 
decay constants and the CKM matrix elements. In both cases, improvement can 
be anticipated with refined lattice QCD calculations. 

The suppression of the —> /r + /r _ branching fractions is characteristic of 

the SM. In particular, scalar contributions from SM Higgs penguin diagrams are 
negligible due to the small size of the muon mass. However, many BSM models 
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can cause the branching fractions to deviate from their SM values. Models with 


extended Higgs’ sectors, for example, can produce significant enhancements in 
the rates of the decays as the helicity suppression is broken (see, for example, 


Refs. 91-94J). The ratio of branching fractions for B° and B® decays to dimuons 
also provides a stringent test of MFV. 

Prior to data taking at the LHC, no evidence for either decay had been found 
and limits on their branching fractions were still an order of magnitude above 


the SM expectations 95, 961. A series of results from LHC experiments 197 


101 significantly restricted the available phase space for BSM theories, giving 


strong constraints complementary to those from searches for on-shell production 
of new particles. In summer 2013 both CMS and LHCb were both able to report 
evidence for the B® —> /x + /r _ decay at the level of four standard deviations (a) 
using their full Run I datasets [102,103 . The experiments exploit multivariate 
event classifiers to optimise the separation of signal from backgrounds consisting 
of muons from different 6-hadron decays. Backgrounds from 6-hadron decays 
where one or more particles were mistakenly identified as muon or where one or 
more particles from a 6-hadron decay was not reconstructed were estimated using 
samples of simulated events whose performance had been corrected to match that 
of data. 

To obtain the best information from the LHC Run I dataset, a simultaneous 


analysis of the CMS and LHCb datasets was performed, giving 104 


B(B"^ /rV) = (2-8 x 10 “ y and 


(19) 


B{B 0 ^ fj+fT) = (3-9 x 10“ 10 . 

These measurements constitute the first observation of the B® —» /x + /i“ decay at 
more than 6a and the first evidence for the B° —> decay at more than 3a. 
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Figure 3: Likelihood contours for the B !? — > /r + /i~ and B° —» n + fi~ branching 
fractions from a fit to the combined LHCb and CMS datasets. The cross indicates 
the best fit point in the two-dimensional plane of the branching fractions, where 
the SM and MFV predictions are also shown. The one dimensional projections 
of the likelihood are also shown. Modified from Ref. |l04 . 

The results, shown in Fig. [3j are compatible with the SM at the level of around 
2 a and put strong constraints on possible BSM scalar and pseudoscalar operators 
as discussed in Sec. [H 

The observation of the B® —> /r + /r~ decay and the first evidence for the B° —> 
/x _ decay represent the culmination of an experimental search that lasted 
three decades. The progress for the other dilepton modes has been somewhat 
less dramatic — as summarised in Table [2j the branching fraction limits are at 
least five orders of magnitude above the SM expectations. Prospects for improved 
measurements are discussed in Sec. [5j 

3.2 Radiative decays 

As mentioned in Sec. [lj measurements of the B — > A s y branching fraction, per¬ 
formed by BaBar and Belle [5j|7| are consistent with the SM expectation [8j. 
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Table 2: Theoretical predictions [84| and experimental results for time-integrated 
branching fractions of B meson decays to dilepton final states, B(B9s —> £ + i~). 
Upper limits are at 90% confidence level. 


Decay 


B°- 

B°- 


e 1 e 


e 1 e 


Prediction 

(2.48 ±0.21) x 10 -15 
(8.54 ±0.55) x 10 -14 


Measurement 

< 8.3 x 10 -8 

< 2.8 x 10 -7 


105 


105 


(1.06 ± 0.09) x IQ’ 1 * * * * * * * * 10 (3.9 + }j) X 10" 10 


B/i+,u~ (3.65 ± 0.23) x 10~ 9 (2.8 ±g;|) x 10 


1-9 


B°^t + t~ (2.22 ±0.19) x 10~ 8 

B° s -A t + t~ (7.73 ± 0.49) x 10~ 7 


< 4.1 x 10 


-3 


104 


104 


106 


No result^ 


1 A limit B(B °—> r + r ) < 5.0% has been derived 107 from ALEPH data 108 . 


Therefore, the main focus in this area has switched to ratios of processes medi¬ 

ated by b — > d'y and b — > sy transitions, CP and isospin asymmetry measurements, 

and measurements of the polarisation of the emitted photon. 

The ratio of branching fractions for b —> d'y and b —> s-y mediated decays 
determines the ratio of CKM matrix elements \Vtd/Vt s \ 2 , and is therefore of in¬ 
terest to search for non-MFV BSM signatures. The measurement of the inclusive 

B(B —>• X c f-j) is very challenging for any experiment, but has nonetheless been 
performed by BaBar 109j. Perhaps more promising is the possibility to use 


exclusive decays, such as B(B° p°'y)/B(B° -A K*° 7) 68,69 . Results have 
been presented by both BaBar iig and Belle nil], and give a precision on 

|Vf,rf/Vts| of about 10%. Further reduction in the experimental uncertainty can 
be anticipated with results from LHCb, which has demonstrated its potential to 


reconstruct B° —> K *°7 112 by making the most precise determination of its CP 


asymmetry. Indeed, CP asymmetries of both inclusive and exclusive radiative 
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Table 3: Measurements of CP and isospin asymmetries in b—> s'f transitions 119]. 
The value of Acp{B —> K* 7) is the average for the B° decay which is much 
more precise than the value for the B + decay. The value of Acp(B —»• X s j) is 



B^X S 7 

B -> K*~f 

Acp 

0.015 ±0.020 

-0.002 ±0.015 

A\ 

-0.01 ± 0.06 

0.012 ±0.051 


b hadron decays offer powerful null tests of the SM 113 , as do isospin asymme¬ 


tries (i.e. differences between charged and neutral B meson decay rates) H 
All such measurements to date are consistent with the SM as shown in Table |3l 

The branching fractions of inclusive and exclusive 6 —> s'y decays are propor¬ 
tional at leading order to the photon dipole operator squared, \C 7 \ 2 ± || 2 , and 
are not sensitive to the handedness of the emitted photon. In the SM, photons 
produced in radiative b hadron decays are almost entirely left-handed due to the 
chiral nature of the charged current interaction. The right-handed component 
is suppressed by the ratio C 7 /C 7 ~ m s /m\ } , with the exact level of suppression 
being mode dependent due to QCD effects |69,116 . In many models that ex¬ 
tend the SM, the virtual particles responsible for mediating the decay have no 
preferred left- or right-handed coupling and the photon can be produced with a 
significantly lower degree of polarisation. Well-known examples include super- 
symmetric models beyond MFV, left-right symmetric models, leptoquarks and 
models with additional gauge-bosons. 

Several methods to measure the photon polarisation have been proposed. One 
of the most promising exploits the interference between B° and B° decays when 
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the hadronic system in the B—tfj decay is accessible to both. In such a case, 


the decay time dependent asymmetry can be written 


A C p[B -> /7](t) = Sj 7 sin(Ami) — Cf 1 cos(Amf) , 


( 20 ) 


where Am is the mass difference in the B°-B° system. The CP asymmetry in 
decay, C/ 7 , has the same sensitivity as that measured with decay-time integrated 
methods, but the coefficient of the sinusoidal oscillation can be written 1117,118] 


S fi = Xf sin(2 ip) sin(2/3) , 


( 21 ) 


where X/ is the C eigenvalue of the hadronic system, tanif> gives the magnitude 


of the ratio of right- and left-handed amplitudes and /3 = arg 


v td v t * b 


is the 


angle of the CKM Unitarity Triangle, which is measured to be sin(2/3) = 0.682 ± 
assuming no NP in B°-B° oscillations. Measurements of Sf~ can 


0.019 


19 


therefore be interpreted in terms of sin(2^), and thus in terms of C^/Cy. 

The coefficient Sk *7 has been measured by BaBar and Belle using B° —)• I \*°7 


decays where K*°—> K^tt 0 , giving 19,119,120 


S K * 7 = -0.16 ±0.22 


( 22 ) 


which is consistent with the SM prediction of —0.02 69 . Measurements of similar 


coefficients in different final states are somewhat less precise, though the results 
for B° —> K®p °7 are competitive 


121 


122 


. Significant improvement in the sen¬ 


sitivity to right-handed currents is a key goal of current and future experiments. 

The B° -> KqTt 0 7 decay is highly challenging to reconstruct at the LHC, but 
the production of all b hadron species opens alternative possibilities. Sugges¬ 
tions to measure the photon polarisation using A® —> 7 decays 


123 


124 


at 


the LHC have proved experimentally difficult due to the small polarisation of 
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. However, the decay appears 


v vx 

ca cb suppresses the 


attractive. Although the small SM value of f3 s = — arg 
Sfj coefficient, the non-zero value of the width difference Ar s in the B®~B® sys¬ 
tem results in sensitivity to the photon polarisation through the effective lifetime, 


112 


or equivalently the A Ar parameter, of B® —> </>7 decays 126 . An experimental 
advantage of such an analysis is that it does not require flavour tagging. LHCb 
has previously shown that it can reconstruct large yields of B J? —> <j)^ decays 
and its first results on the effective lifetime are keenly anticipated. 

Another way to probe C' 7 jC 7 is through the photon direction with respect to 
the plane defined by the 7r + 7r - system in B + —> decays [i27jjl29j 


The so-called up-down asymmetry of the photon with respect to this plane is 
proportional to the photon polarisation. The constant of proportionality suffers 
large hadronic uncertainties, that can, however, be controlled to some extent from 
data on B + — > J/ip K + ir + ir~ decays and other information concerning the Kirn 
system. LHCb has performed a first measurement of this up-down asymmetry 
using its full Run I dataset |130], As shown in Fig. [4j the data are split into 
four regions of A' + 7 t - 7t + mass defined by the known A + 7 t _ 7t + 

resonances. Whilst non-zero photon polarisation is observed at the level of 5.2cr, 
when combining the four regions of M(K + tt~'k + ), a deeper understanding of 
the structure of the A + 7r _ 7r + system is needed to determine with this 

approach. 

The B° —> I\* 0 £ + £- decays can also be used to determine C' 7 /C 7 at low dilepton 
invariant masses, where virtual photon contributions are expected to dominate. 
Study of the angular distributions allows the determination of the parameter 


Ay defined in Eq. (12). At the limit q 2 —> 0 this quantity is directly sensitive to 
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Figure 4: Top: background subtracted K + tt~tt + invariant mass distribution 
in B + —> /\' + 7r _ 7r + 7 decays. Bottom: up-down asymmetry A U( j in bins of 
Modified from Ref. 


130 


C' 7 /C 7 (considering complex valued Wilson coefficients, another observable 
probes the relative phase between C' 7 and C 7 ). LHCb has studied B° —> I\*°e + e~ 


decays in the low q 2 region 131,132 , obtaining 


= -0.23 ±0.23 ±0.05 


( 23 ) 


for 0.002 < q 2 < 1.120 GeV 2 /c 4 from the Run I dataset. 
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3.3 Semileptonic b -h s£+r decays 


The LHC data has led to a wealth of results in semileptonic b —> s£ + l~ decays. In 
the following subsections results on branching fractions and rate asymmetries are 
first discussed, followed by considerations on analyses of angular distributions. 
Branching fractions of semileptonic b -a s£ + £~ decays In contrast 
to the case of b — > sj decays, most of the experimental work on the b —> s£ + £~ 


process has been with exclusive final states. Results on the inclusive decay 133 


134 do not yet reach high precision, but will be an important topic in the future. 


The LHC data have, however, led to a large increase in the yields of certain 
b —> s£ + £~ decays, in particular those with a dimuon pair in the final state. This 
has led to increasingly precise determinations of the branching fractions of the 
B—> B—> and B® —> <;A/z + /i _ decays 

semileptonic b baryon decay A® —> A^i + 


135 


138 


, as well as of the 


139 


140 


. The theory input for the 


baryon decays 141,142|, in particular knowledge of the form-factors, is less well 
advanced than for the mesons, so they are not discussed further in this review. 

The experimental results on the heavy-to-light branching fractions are now 
much more precise than the corresponding theoretical predictions, with further 
improved measurements anticipated in the coming years. This is illustrated in 
Table[4j where measurements of B(B + —> K + fr + fi~) are compared to theory pre¬ 
dictions in the low and high q 2 regions. SM predictions for the branching fractions 
are sensitive to hadronic uncertainties in the form factors, which typically lead to 
uncertainties of 0(30%) on the SM predictions. Progress from lattice QCD has 
improved the precision in the high q 2 region above the i/^(3770) resonance, but 
the sensitivity to BSM physics still remains limited by the uncertainty on the SM 


predictions. 
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Table 4: Branching fraction of the B + —> K + fi + /i decay in selected q 2 bins, 

fe 2 


mmiTmax] (in GeV 2 /c 4 ), from LHCb 143 compared to SM predictions us¬ 
ing light-cone-sum-rule (LCSR) 138] and lattice |40j calculations of form factors. 
The first uncertainty on the experimental results is statistical and the second is 
systematic. 



LHCb 

SM (LCSR) 

SM (lattice) 

[1.6] 

[16,18] 

(1.21 ±0.09 ±0.07) x 10“ 7 

(0.35 ±0.04 ±0.02) x 10“ 7 

(1.75™) xl0“ 7 

(0.33™) xlO- 7 

(1.81 ±0.61) x 10“ 7 

(0.39 ±0.04) x 10“ 7 


In order to increase sensitivity to BSM physics, it is useful to study observables 
in which the effects of form factor uncertainties are reduced. Two such quantities 
are the CP asymmetry between B and B decays, Acp, and the isospin asymmetry, 
A\ between B + and B° decays. In the SM, the CP asymmetries of B—> 
and B -> Kfr + fi~ decays are tiny, 0(1O~ 3 ), due to the small numerical size of the 
product of CKM elements V ub V* s compared to V tb V t *. Extensions of the SM can 
provide new sources of CP violation, and thus Acp constitutes a null test of the 
SM, where any visible direct CP violation would be evidence for BSM physics. 
The latest results, all consistent with the SM expectation of close to zero, are 
summarised in Tabled 

The decays appearing in the isospin asymmetry only differ by the flavour of 
the spectator quark in the B meson (u quark for the B + meson and d quark for 
the B°). In the framework of the effective Hamiltonian, Mi differs from zero due 
to isospin breaking effects in the form factors, from annihilation and exchange 
amplitudes, and from spectator scattering, where a virtual photon is emitted. 
The isospin asymmetry is expected to be ~ —1% at large q 2 in the SM 55 and 
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Table 5: Measurements of CP and isospin asymmetries in B —> 
transitions at low and high values of q 2 135, 145 . Units of GeV 2 /c 4 for q 2 


are implied. Only results from LHCb are included since earlier, less precise, 
measurements used a different q 2 binning scheme. 

B^Kh+ij,~ B—> K*n + ^~ 


q range 


1 . 1 - 6.0 


15 . 0 - 22.0 


1 . 1 - 6.0 


15 . 0 - 19.0 


Acp 

Ai 


0.004 ± 0.028 


- 0.005 ± 0.030 


- 0 . 10 +°; g ®± 0.02 - 0.09 ± 0.08 ± 0.02 


- 0.094 ± 0.047 - 0.074 ± 0.044 


0.00 ± 0.02 0 . 06 ± 0.02 


for B —» I\* decays to grow to « ±10% as q 2 tends to zero 114,144 . Using 


its full Run I dataset, LHCb found A\ in B -± K and B —> K* decays to be 
consistent with zero across the full q 2 window 1 135] (see Table [ 5 ]). 

The ratio of rates for b —> d£ + £~ and b —> s£ + £~ processes is sensitive to 
\Vtd/Vts\ 2 - The observation by LHCb of the B + —> 7 r + /r + /u decay gives 1146 


B(B + —> 7 T + qi + fi ) = (2.3 ± 0.6 (stat) ±0.1 (syst)) x 10 8 , 


(24) 


\V td /V t s\ = 0.266 ± 0.035 (stat) ± 0.003 (syst) , 
where the uncertainty on | V t d/Vt s \ due to knowledge of the form factors, estimated 
to be 5.1%, is not included in the result. Further improvements, and observations 
of more b — > d£ + l~ decay modes, are anticipated in the coming years. 

Angular ANALYSES of b —> s£ + £~ DECAYS An angular analysis of B —> 
K(. + H decays provides a simple null test of the SM |147 . The angular distribu¬ 
tion can be described by a single angle 9g 

— - : —— — = t(1 — -^h)( 1 — cos2 ® l ) + + ^fb cos Qt , (25) 

1 d cos 4 2 


with a constant term, T H /2, and a forward-backward asymmetry, A FB , linear in 
cos Both and A FB are small within the SM, for l = e, /r, and therefore can 
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signal the presence of BSM physics. In particular these terms are sensitive to con¬ 


tributions from new scalar, pseudoscalar and tensor operators. LHCb has made 
precise measurements of these parameters in the decay B + — > J\ + /r + / u~ 1 148 . 
The measurements are consistent with H FB = 0 and o, expected in the SM. 

The angular distribution of the B° —> A'*°/x + /i _ decay, with A"* 0 —> K + n ~, 
is more complicated, and can be described by three angles: 9p, which is defined 
by the direction of the (/i _ ) with respect to the B° ( B °) in the dimuon rest 
frame; 9k, which is defined by the direction of the kaon with respect to the B° 
(B°) in the K*° (K*°) rest frame; and </>, the angle between the plane containing 
the n + and and the plane containing the kaon and pion. The differential 
decay rates in terms of these angles and the dimuon invariant mass squared, for 
B° and B° decays, are given by 


d 4 T[A° —► K*°n + fi-} _ 9 


d cos 61 d cos 9k d</> d q 2 327T z - 

d 4 T[A° A* W] 9 


yj^Mcose^cos Q k ,4>) 


(26) 


y Ji(q 2 )fi(cos 6i, cos 9k, 4>) ■ 


d cos 9 1 d cos 9k d</> dq 2 32ir 

Here, the /*(cos 9g, cos 9k, <t>) originate from spherical harmonics and the Jj and J, 
are bilinear combinations of K*° decay amplitudes and Aq’ R ) 

The CP averaged observables, 


149 


Si — (Ji + Ji) 


'd(r + r) 

dq 2 


(27) 


depend on the underlying short distance contributions from C 7 ± C' 7 , C 9 ± C' 9 
and C 10 ± C [ 0 , for CP odd (even) components. They can be related to the 
observables discussed in Sec. 


2.3 


for example by S 3 = (1 — F L )A^/2 , S 4 5 = 


y/F L (l — FjJPi 5- In addition, CP violating observables (oc Ji — Ji) can be ob¬ 
tained from the angular distributions, including several with high BSM sensitiv¬ 
ity [l49j- 152 . 
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Theory ^■Binned 
-LHCb ATLAS -«-CMS 



Theory ^■Binned 
-LHCb ATLAS -"-CMS 




HH= 



15 20 

q 2 [GeV^c 4 ] 


15 20 

C 7 2 [GeV 2 /c 4 ] 


Figure 5: (Left) longitudinal polarisation fraction, TL, of iL*° mesons produced 
in P°—^ decays. (Right) dimuon system forward-backward asymmetry, 

Ape- Results from ATLAS [l53|, CMS [137 and LHCb |l36j are included. The 


data are overlaid with an SM prediction 154 


Prior to data taking at the LHC, with the relatively modest samples of B° —> 
decays that were available, it was not possible to determine all of the 


terms of Eq. (26). Instead, partial angular analyses of the decay were performed 


using single angle projections. ATLAS, CMS and LHCb have also performed 

similar analyses 136,137,153 giving sensitivity to Fl, Afb (see Fig. [5]) and, in 
( 2 ) 

the case of LHCb, . Each of these observables is consistent with the SM 
expectation. 

LHCb has also measured two of the “optimised” observables discussed in 

In the low q 2 region, there is a large local discrep- 


Sec. 


2.3 


P \ and Ei 


155 


ancy between the data for and the SM expectation at the level of 3.7cr, as 
shown in Fig. [6j This is discussed further in Sec. [4} 

It is expected that full angular analyses of the B° —> K*°qi + decay should 
be possible with the full Run I datasets of the LHC experiments. As the analyses 
become more precise, and more complex, it will also be important to account for 
the contribution from Kn S-wave under the K*° peak 1157 -1591. It will also be 
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Figure 6: Observables P' x and PL measured by LHCb 


155 


in B u -> K*°n + n 


decays. The data are overlaid with an SM prediction 156 


possible to determine CP asymmetries for each of the angular terms; these can 
then be used to constrain the imaginary parts of the Wilson coefficients. Other 
B—> V£ + £~ decays, such as B® —> (fru+fj,-, will provide additional constraints. 


3.4 Non-universal lepton couplings 


In the SM, with the notable exception of the Higgs boson, particles couple equally 
to the different flavours of lepton. The ratio of decay rates 


Rh = r[£?-> ify+/z- 


(28) 


T[B^ He+e -] 

where H = K, K *, X s , etc., is therefore expected to differ from unity only due 
to tiny Higgs penguin contributions and phase space differences 1601. Using the 


full Run I dataset, LHCb measured in the q 2 range 1 < q 2 < 6GeV 2 /c 4 161 


R k [ 1, 6] = 0.745 (stat) ± 0.036 (syst), 


(29) 


which differs from the SM expectation of Rk = 1.0003 ± 0.0001 |147 by 2.6cr. 
Although not yet at the level of significance that qualifies as “evidence”, this 
result has prompted theoretical speculation concerning possible sources of lepton 
non-universality, as discussed further in Sec. [|J Since results from BaBar on 
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also hint at violation of universality, this is a highly 


3.5 Null tests 

Null tests of the SM, i.e. searches for signals that are absent or vanishingly small 
in the SM, are valuable for several reasons. Observation of such a process would 
not only provide a smoking gun signature of BSM physics, but would also indicate 
in what way the NP should be accommodated within the operator basis discussed 


in Sec. 2.1 In the absence of signals, limits can be placed on the contributions 


of the additional operators which, when sufficiently stringent, justify the use of 
a restricted set of operators in model-independent analyses. 

The tests of lepton universality discussed above fall into this category of null 
tests. It is also important to explore the possibility of lepton flavour violation 
(LFV), and lepton number violation (LNV) in rare b hadron decays. The ob¬ 
servation of neutrino oscillation demonstrates that lepton flavour is not an exact 
symmetry of nature, but if the SM is minimally expanded to allow neutrino mass, 
rates of processes with charged LFV remain unobservably small. Charged LFV, 
and LNV, does however arise in many BSM theories. While strong limits exist 
from searches for rare muon and tau decays (for LFV) |164| and for neutrinoless 
double beta decay (for LNV) [165 , there are models which respect those limits 
but nonetheless produce observable signatures in b hadron decays. 

One of the most powerful of such LFV searches is for the —y decay. 


The experimental limits have been improved by LHCb to the level of < 10 8 166 


The limit on the B° —»• decay is four orders of magnitude more stringent 


than those on B° —> and B° —> /r ± r T 167 . For semileptonic decays, 
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the limits on the B + —> -k + e±/jX and B + —> K + e ± n T decays, at the level of 

10- 7 168,169 , are similarly much more stringent than those on decays involving 


r leptons |170|. Since the operators can be written as being independent for each 
pair of leptons, it is important to improve all limits. These LFV semileptonic 
decay modes have not yet been investigated at the LHC. 

Several searches have been performed for LNV in b hadron decays into final 
states containing a pair of same sign leptons, B + —> M~l + £' + . The strongest 


limits, from LHCb on B(B + — > 7 r _ /r + /i + ), are at the level of 10 -9 171 . These 
are complemented by limits on numerous modes with different hadronic systems 
= 7 t~, p~, K~,K*~ , D~, D*~ , Dj ,... and in which one or both leptons may 


be an electron, all at the level of 10 6 or below 172-176 . Limits also exist at the 


10 6 level 177 on several b hadron decays that violate both baryon and lepton 
number, but these have not yet been explored at the LHC. 


The most recent LHCb analysis of B + —> 7r“^ + /i + 171 sets limits on the 
branching fraction as functions of the mass and decay time of the 7 pair. 
These results are of interest to probe models where the decay is mediated by the 
on-shell production of a Majorana neutrino, which could be long-lived. Similar 
experimental techniques can be exploited to search for long-lived particles (A") in 
B—t KX with X—> l + i~ and similar decays. Such signatures are predicted in a 
range of theories, that are generically referred to as “dark sector” models 1178 


Although these decays have not yet been investigated at the LHC, the possibility 
to search for particles that travel 0(1 m) before decaying makes these probes 


complementary to other searches for new light resonances. 
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4 Interpretation 

4.1 Wilson coefficient fits 


The observation of the B° s —> y + y decay, described in Sec. 3.1, puts strong 


constraints on scalar and pseudo-scalar operators (O^ and O^}). In BSM models, 
the Bg —> £ + £~ branching fraction is enhanced or suppressed by the ratio 
B(E° a ->£+£-) 


0+p , = |1 —0.24(Cfo F —Cfo)—y£(Cp —Cp)f+ |y£(Cg —Cg )f , (30) 

where y^ = 7.7 and y e = {m ll /m e )y pi = 1.6 x 10 3 . At ler, the current experimental 
measurements imply [571 


|Cp — Cp\ < 0.3 and ^-^<0.1. 


(31) 


The constraints for dielectron decays are somewhat weaker, 


I Cg — Cg [ + |Cp — C'p \ z !.3 , 


(32) 


and substantial room still exists for ditau decays. Barring fortuitous cancella¬ 
tions that can happen if C s = C' s and C P = C'p, visible effects from operators 
O g and Op to dimuon and dielectron decays can be neglected. Large acciden¬ 
tal cancellations are also excluded by B —> K£ + £~ decays, which constrain the 
combinations |C S + Cg| and |C P + Cp|. Contributions from tensor operators, are 
also constrained by B —> K£ + £~ decays. In particular, the small size of the F p 


term in the B + —> K + y + y angular distribution 33 , leads to a bound of 


|Ct| 2 + |Ct5| 2 < 0.5 . 


(33) 


Therefore, it is appropriate to focus on the operator basis of Eq. ([2]), and consider 
what the current measurements tell us about BSM contributions. 

The rate of B —> A s y decays is consistent with the SM and places constraints 
on the size of BSM contributions to the electromagnetic dipole operators Cy\ 
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The LHCb measurement of the up-down asymmetry in B + —> Ji + 7 r _ 7 r + 7 decays 


shows that the emitted photons are polarised, but further work is needed to 
interpret the asymmetry. Moreover, measurements of Sk *7 from BaBar and 
Belle are consistent with the SM expectation that C 7 3> C 7 , but are not yet 
sufficient to rule out a sizable right-handed polarisation. Figure [7] illustrates 
the current experimental constraints on C 7 after fixing all of the other Wilson 
coefficients to their SM values. A future precise measurement of at low q 2 
in the B° —> K*°e + e~ decay in combination with reduction on the uncertainty 
on Sk* 7 would rule out a large right-handed contribution 1179 . 

While C 7 and C 7 are consistent with their SM expectations, the situation with 
c!j \ q is more interesting. Measurements show a trend for the rates of the B —> 
K (*) fr + fr and B® —> (j)fi + n decays to be below their SM expectations at both 
low and high q 2 . The angular observable P§ at low q 2 also appears to be discrepant 
with respect to the SM, although the other angular observables in the B° —> 
decay are reasonably consistent with their SM expectations. Figure [ 8 ] 
illustrates how these measurements relate to BSM contributions to C 9 and C 9 . 
In general the data, while still consistent with the SM, are best described by a 
destructive BSM contribution to C 9 , which both reduces the branching fraction 
of the semileptonic b —> sfi + decays and also modifies the angular distribution 
of the B 0 -^- decay to be more consistent with the observed value of 

at low q 2 . Global fits to b —> s data favour Cg^ ~ — 1 with other NP parameters 
consistent with zero or small additional contributions to C 9 or C 10 
Similar conclusions have been obtained from analyses that differ in statistical 
treatment, theoretical input (form factors) and in the treatment of systematic 
uncertainties (incorporation of power corrections). A reduction in the rate of 
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b —> s[i + [r~ processes could also be achieved by a contribution to C w with the 

opposite sign to the SM (C^, p > 0), as shown in Fig. [ 9 } However, large BSM 
contributions to C 10 are disfavoured by the measured branching fraction of the 
B° s —> fr + fi decay. 

Current fits exhibit very little sensitivity to the phases of the Wilson coefficients 
This situation could be improved through the measurement 


Cf and C'fJ 


154 


of CP asymmetries in the B° —> K* 0 /j, + n~ and B® —> (j)fr + fr~ angular distri¬ 
butions. Some of these asymmetries are naive T-odd and hence do not require 
sizable strong phase differences between the decay amplitudes in order to be 
non-vanishing. 


The measurement of Rk, Eq. (29), raises the question of whether to consider 


only universal contributions to C 9 \ 0 . Considering flavour-dependent Wilson co¬ 
efficients, the data imply 


0.7 < Re[X e - <1.5, 

where A £ = C 9 NP£ + C' 9 e - (Cf 0 P£ + Cf 0 ), £ = e,fi. 


(34) 


The anomaly could be caused by BSM physics either suppressing the dimuon 
mode, enhancing the (currently less constrained) dielectron channel, or both. 
Interestingly, whilst the branching fraction of the B + —> /i + /i + /z - decay is below 
the SM expectation, that of the B + —> K + e + e~ decay is consistent with the 
prediction. The ratios Rk and Rk* are interesting null tests of the SM because 
their theoretical predictions are free of hadronic uncertainties. A measurement 
of Rk*, or of the ratio of inclusive b s£ + l~ branching fractions, Rx s , in 
combination with the measurement of Rk could be used to separate non-universal 


BSM contributions to C g 10 from those to C 9 10 184 . For semileptonic decays 


with a ditau, or dineutrino, pair, the current data leave ample room for BSM 
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Figure 7: Dependence of the observables B(B —> X s j), Sk*j and in the q 2 
range 0.002 < q 2 < 1.120 GeV 2 /c 4 on Re(C 7 ) and Irr^C^), computed using the 
EOS flavour tool 1621 with all other Wilson coefficients fixed to their SM expec¬ 
tations. The contours (red) indicate the 68% and 90% intervals on experimental 
determinations of the observables. The lines indicate the theoretical uncertainty 
on each observable. 

contributions. In most models, however, correlations exist between the different 
lepton flavours and indirect bounds on the couplings can be derived from the 
dimuon and dielectron data. 

Global analyses of the b —> s data have also been explored using the SU(2 )l 
invariant basis described in Sec. |2.1| 1 183,185|[l86j. These global fits have pref¬ 
erence for a BSM contribution to Cll = Cq — C io, with left-handed muons, as 
opposed to C*lr = Cg + Cio- 

Similar analyses of results on b —> d decays are also possible, though the preci¬ 
sion of the current data is not sufficient to obtain useful constraints. Comparison 
of b —> d with b —> s data allows to search for non-MFV BSM signatures. Since 
violations of lepton universality may themselves be non-MFV signatures, if the 
hints for non-universality discussed above persist, it will become even more im¬ 
portant to improve the precision of the constraints in the b —>• d sector. 
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^ FB I 1 ’®1 
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Figure 8: Top: dependence of the angular observables A FB , P§ and A T in 
the B° —> K* 0 ir + /r~ decay in the q 2 range 1 < q 2 < 6GeV 2 /c 4 on Cg F and 
C 9 - bottom: dependence of the branching fractions B(B° —> K*°fi + g~) and 
B(B + -» K + g + g~) in the q 2 range 1 < q 2 < 6GeV 2 /c 4 ; all computed using 


the EOS flavour tool 62 with all other Wilson coefficients fixed to their SM 


expectations. The contours indicate the 68% and 90% intervals on experimental 
determinations of the observable. The lines indicate the theoretical uncertainty 
on each observable. 

4.2 Limits on NP scales 

In the SM, FCNC amplitudes for |A.B| = |A5| = 1 transitions are suppressed 
because they occur through loop effects involving the weak scale (g 2 /(47r) 2 and 
1/M^r), and also due to the smallness of the relevant CKM matrix elements, as 
seen in Eq. 0. As discussed in Sec. |2.4[ NP models may share some or all of 


these features. The assumed amount of suppression of BSM amplitudes impacts 
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Figure 9: Top: dependence of the angular observables A FB , P 5 and A T in 
the B° — y K*°/x + ^~ decay in the q 2 range 1 < q 2 < 6 GeV 2 /c 4 on C^q 3 and 
C[ 0 ; bottom: dependence of the branching fractions B(B° —> K*°fr + fi~) and 
B(B + — » K + /i + fi~) in the q 2 range 1 < q 2 < 6 GeV 2 /c 4 and of B{B° s ^f /x + ^ _ ); all 
computed using the EOS flavour tool 1621 with all other Wilson coefficients fixed 
to their SM expectations. The contours indicate the 68 % and 90% intervals on 
experimental determinations of the observables. The lines indicate the theoretical 
uncertainty on each observable. 

on the limits on the corresponding scales. 

The best fit values for C^f 0 can be interpreted in terms of a BSM scale Anp, 


Anp x ^/|C^f 0 | 




4 ttV 2 m w _ gg rpev (generic tree level), 
- 2Tev 

\/2M\y/e = 400 GeV (MFV, weak loop). 


(35) 


Thus, the value of |Cg | = 1 obtained in Sec. 4.1 corresponds to Anp ranging 
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from ~ 400 GeV to ~ 36 TeV depending on the model-dependent suppression of 


the BSM amplitude. For |AF>| = |AD| = 1 transitions, the CKM suppression in 
the SM is stronger and an analogous bound would imply a stronger constraint 
on Anp in models without flavour-suppression by \f\Vts/V t d\ ~ 2. 

If explicit flavour factors, AA*, are introduced for the BSM contribution, rare 
decays provide constraints on the combination AA*/A^ p . In contrast, B~B mix¬ 
ing constrains the combination (AA*) 2 /A^ p . Similarly in the SM, the loop con¬ 
tribution to mixing is proportional to Gf| V^V^* | 2 . Due to the stronger CKM 
suppression and in view of the strong constraints on CP violation in mixing, 
mixing bounds on Anp are typically more powerful than those from rare decays 
with two important exceptions. The first is where the A B = 1 BSM amplitude 
arises at tree level but the A B = 2 BSM amplitude is loop induced, as in lepto- 
quark models. Second, if the flavour suppression in the BSM model is sufficiently 
strong, i.e. if AA* is small, rare decays provide more stringent constraints on Anp. 
It should also be emphasised that due to the different dependence on scales and 
couplings, combinations of measurements of A B = 1 and A B = 2 processes can 
be used to fix both dimensionful and dimensionless BSM parameters. 

4.3 Impact on model building 

The available data on rare b hadron decays put strong constraints on BSM contri¬ 
butions to the amplitudes, and this has significant impact on the range of models 
that can be considered. As discussed above, BSM effects in the Wilson coeffi¬ 
cients of the semileptonic vector (C 9 ) and axial-vector operators ((7 10 ) are now 
limited to be no larger than roughly a third of their respective SM values, with 
constraints on the chirality-flipped coefficients (Cg 10 ) of similar size. The Wilson 
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coefficients of non-SM operators (scalars, tensors) are also strongly constrained, 


as shown in Eqs. (31) and (33). Note that these limits are for muons, and there 


is room left for possible signals in decays into other lepton species. 

A common source of (pseudo-)scalar operators in BSM models is Higgs-induced 
penguins. These processes are Yukawa-dependent and hence couple differently to 
the different flavours of lepton. In the minimal supersymmetric standard model 
(MSSM) with MFV, the introduction of a second Higgs doublet can lead to an 


enhancement of Cg and Cp 


93 


94 


Cg ~ —Cp oc me tav? f3/m 2 A . 


(36) 


Here, tan /3 is the ratio of the vacuum expectation values for the two Higgs dou¬ 
blets present in the theory and rriA is the mass of the CP-odd pseudoscalar Higgs. 
A sizable value of tan/3, of around fifteen, can overcome the small muon mass 
factor. The constraints on Cg and Cp arising from the branching fraction of 
the B® —> /r + /Y decay are therefore able to rule out significant amounts of the 


phase-space of MSSM models 187 


While the overall picture is one of consistency with the SM, there are some 
hints of anomalies in the b —> s data, which if taken at face value exhibit quite 
interesting features: a) a preference for sizable NP in leptonic vector couplings, 
with axial-vector NP not larger in magnitude; b) a preference for lepton non¬ 
universality; c) the possibility of right-handed currents. If future measurements 
substantiate any of these, there will be huge implications for model building. To 
accommodate such features requires models that go beyond the most common and 
simple solutions to the hierarchy problem. For instance, b) and c) directly imply 


a non-MFV flavour sector, while b) also suggests lepton-flavour violation 80 


It should be noted that the significant level of non-universality hinted at in the 
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data, if confirmed, would rule out many proposed SM extensions including the 


MSSM with R-parity conservation |188 


Feature a) has inspired model-building with Z'-extensions to the SM. As ar¬ 


gued in Sec. 2.4, many models with possibilities for large Z-penguins predict 
the hierarchy between NP in axial-vector and vector couplings to be the other 
way around. However, one popular approach that prefers large vector NP is the 
SU(3) x SU(3) x U(l), or “3-3-1”, model 1189,1901, which can also accommodate 
feature c). A model which can explain all of a), b) and c) is the gauged L T ~L^ 


extension of the SM |76,77 . A variant of the latter with an additional Higgs- 


doublet 191 can also explain the 2.5<r hint from CMS of Higgs boson decay to 




192 . In this model there are constraints on the ratio of the mass of the 


Z' and of its coupling, mz' / g ', in the few TeV range. 

Leptoquark models offer a natural framework to accommodate lepton non¬ 
universality. Choosing an SU(2)L-triplet leptoquark that couples to muon dou¬ 
blets one obtains a model that induces at tree level [33,! 193 


pNP(i _ 7T V2 


/-iNP/i _ /inrji _ 

°9 — °10 — 


(37) 


a e V fb V* 2M 2 Gf ’ 

where M is the mass of the leptoquark. With benchmark values to explain the 
measurement of Rk given in Eq. (29), M 2 ~ (48 TeV) 2 . Viable flavour 

structures for the leptoquark couplings X q £ can arise in models with partial com¬ 


positeness 193 . Dineutrino operators are induced, which enhance the branching 
ratios of B —> and B —>• decays by a few percent. Corrections of a 

few percent to C 7 can also arise. The relation between axial-vector and vector 


coupling of Eq. (37) can be broken if more than one leptoquark is introduced. 
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The leptoquarks may be searched for through their decays 

( i i 2 / 3 — > t v , (/> -1//3 —>• b v ,t /i _ , 0 _4//3 —tb[j,~, (38) 

where the final state particles must be from different generations since leptoquarks 
carry two generational indices. Such distinctive signatures will however only be 
visible in experiments in the leptoquarks are sufficiently light. 

5 Summary and outlook 

Run I of the LHC has led to a substantial improvement in precision in several 
key observables among rare decays of b hadrons. Particularly notable are the 
first observation of the very rare —>• decay and the wide range of kine¬ 

matic observables now studied in B° —>• iL*°/r + /r _ decays. The results remain 
broadly consistent with the SM, but deviations are present at an intriguing level 
of significance. In the light of this situation, and the high sensitivity to BSM 
physics provided by rare b hadron decays, it is essential to continue to improve 
the precision. 

This presents challenges for both theory and experiment. One important task 
for theory is to reduce uncertainties that arise due to hadronic effects in the 
decays. This includes improving precision in decay constants and form factors, 
for example through refined lattice QCD calculations. In the case of b—> s£ + £~ 
decays, other important issues are to understand power corrections and resonance 
contributions. To address these issues will require improved computations as well 
as detailed studies of data on specific observables and fits. The data will allow not 
only to determine certain hadronic parameters, but will also provide consistency 
checks and allow theoretical methods to be refined. 

On the experimental side, the situation varies between b hadron decay modes. 
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As shown in Table | 2 j the precision of branching fraction measurements of —> 

£ + £~ decays is not yet at the level of the theory predictions; only for B® —y p + 
is it within a factor of 3. Improvements in these measurements can be achieved 
with a combination of experimental facilities. Data taking at the LHC, in Run II 
and beyond, will provide large increases in the yields of dimuon decays. The 
experiments will benefit both from the increased cross-section for production of 
b quarks, which is expected to scale approximately linearly with the centre-of- 
mass energy, and from the high luminosity. The increased yields will allow not 
only improved branching fraction measurements, and a determination of their 


ratio, but also a study of the effective lifetime in B ^ decays 194 . The 

LHC experiments may also be able to improve on the existing limits for the 
dielectron and ditau modes, but substantial improvement is more likely to be 
possible at the Belle II experiment 195 and at a very high luminosity e + e _ —> Z 


factory 196,197 , which have an experimental environment more suitable to these 
modes. 

Radiative b hadron decays provide the potential for significant future improve¬ 
ment in the knowledge of right-handed contributions to the b —> sq dipole am¬ 
plitude. To achieve this it will be necessary to use all of the methods most 
sensitive to the photon polarisation, since they provide complementary informa¬ 
tion 


179 . These methods include time-dependent asymmetries in B° —> K^ir 0 ^ 


and Bg —> 7 decays, up-down asymmetries in B + —> K + TT + iT~'y decays, and 

angular asymmetries in B° —> K*°e + e~ decays. Improved searches for CP 
violation in both inclusive and exclusive processes, as well as tests of MFV 
from B(B° —> p 0, y)/B(B 0 —> K*° 7 ) will also be important. There are excel¬ 
lent prospects for progress in all of these areas at both LHCb, including its 
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upgrade |l98 , and Belle II. 

Data from the LHC have allowed dramatic progress in experimental investi¬ 
gations into semileptonic b—> si + l~ decays. The precision of the measurements 
of the differential branching fractions as a function of q 2 is now good enough 
that the focus is primarily on asymmetries, including CP asymmetries, isospin 
asymmetries, and lepton universality violating differences in rates. With the in¬ 
creasing yields available in decay modes such as B° —> full angular 

analyses are expected to become mandatory. It will also be possible to fit for CP 
violating angular observables. The reconstruction of semileptonic modes with 
dielectron and ditau pairs is more challenging. The immediate objectives are to 
search for violations of lepton universality in a wide range of modes, including 
B° —> I\*°e + e~ and B®-^ <f>e + e~ , and to make first efforts at searching for the di¬ 
tau modes. In addition, high luminosity e + e _ experiments are expected to reach 
interesting sensitivity to the B —> decays, which provide complementary 

information on BSM physics since they are sensitive to other Wilson coefficients. 

The main goal of ongoing investigations into rare b hadron decays is to first 
uncover evidence of BSM physics, and then to deduce its nature. As discussed 
in this review, the quantity and quality of the data being produced by the LHC 
provide exciting prospects. The hints of BSM physics in the data analysed so far 
have led to new directions in model-building with the desired phenomenological 
features. To exploit fully the data that will be available from the LHC and from 
future experiments will require increased effort to identify and interpret patterns 
in the short-distance coefficients. This will provide fantastic opportunities to 
learn about the problems and puzzles in fundamental physics that remain in the 


SM. This challenge can best be confronted by both theory and experiment in 
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collaboration. 


Note added in proof 


After the completion of this review, preliminary results from the full analysis of 
the LHCb Run I data on the angular distributions of the B° —> K*°fi + fj,~ decay 
became available jl99 . Careful study will be required before the detailed im¬ 
plications of these results are understood; however, consistency with the earlier 
results based on a subset of the data [136.1551 is seen and therefore the interpre¬ 
tation of the new data is likely to proceed along similar lines as discussed in this 


review. 
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